Is Z+(4430) a loosely bound molecular state? 



o 
o 

(N 

> 

O 

:^ 

o 
m 



Oh- 
I ■ 

Oh; 



(N 

> 

o 



o 



X 



Xiang LiuiQ Yan-Rui Liu^Ill Wei-Zhen Dengi0 and Shi-Lin Zhu^ 

^Department of Physics, Peking University, Beijing 100871, China 

^Institute of High Energy Physics, P.O. Box 918-4, Beijing 100049, China 

(Dated: February 2, 2008) 

Since ^"'"(4430) lies very close to the threshold of D*Di, we investigate whether Z+(4430) could 
be a loosely bound S-wave state of D*Di or D* D'l with J^ = 0~, 1~,2~, i.e., a molecular state 
arising from the one-pion-exchange potential. The potential from the crossed diagram is much larger 
than that from the diagonal scattering diagram. With various trial wave functions, we notice that 
the attraction from the one pion exchange potential alone is not strong enough to form a bound 
state with realistic pionic coupling constants deduced from the decay widths of Di and D[ . 

PACS numbers: 12.39.Pn, 12.40.Yx, 13. 75. Lb 



I. INTRODUCTION 

Recently Belle Collaboration observed a sharp peak 
in the ir'^ip' invariant mass spectrum in the exclusive 
B -^ Kn'^ip' decays with a statistical significance of 7a 
[l|. This resonance-like structure is named as ^+(4430). 
The fit with a Breit-Wigner form yields its mass m — 
4433 ± 4(stat) ± l(syst) MeV and a narrow width T = 
44+^3 (stat)15°(syst) MeV. It is very interesting to note 
that the width of ^+(4430) is roughly the same as that 
ofDi. 

The product branching fraction is measured to be 
B{B -^ KZ+{AA30)) ■ S(Z+(4430) -^ tt+V'') = (4.1 ± 
l.O(stat) ± 1.3(syst)) x 10"^ 1]. For comparison, we list 
the production rate of X(3872) and F (4260) in B decays. 
From Ref. 2] we have 

B{B- -^ K-X{3872)) ■ B{X{3872) ^ J/^tt+tt") 
= (1.28 ±0.41) X 10"^ 

and from Ref. [31] 

B{B- -^ K-X{3872)) ■ B{X{3872) -^ J/^tt+tt") 
= (10.1 ±2.5 ±1.0) X 10"^. 

For y(4260), Babar Collaboration gave the upper limit 
of the branching fraction Q 

BiB- -^ K-Y{A260)) ■ ^(^(4260) -^ J/^Tr+Tr") 
< 2.9 X 10"^ 

It's plausible that (1) B{B -^ KZ+{U30)) is compa- 
rable to both B{B- -^ K-X{3872)) and B{B- -^ 
ii'~y(4260)); (2) 71+?/)' is one of the main decay modes 
of Z+(4430) if it is a resonance. 

The peak Z+ (4430) inspired several theoretical specu- 
lations of its underlying structure. Rosner suggested that 
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Z+(4430) is a S-wave threshold effect because Z+(4430) 
lies close to the D*(2010)£)i(2420) threshold [1]. The 
production mechanism was speculated as follows. The b 
quark first decays into a strange quark and a pair of cc 
while a pair of light quarks are created from the vacuum. 
In other words, B meson first decays into a K and a pair 
of D mesons. Then the D meson pair re-scatters into 
TT+V'- He also suggested other possible charged states. 

Maiani, Polosa and Riquer identified this signal as the 
first radial excitation of the tetraquark supermultiplet 
to which X(3872) and X(3876) belong Q. With their 
assignment the quantum number of Z+(4430) is J^*^ = 
1^ where the C-parity is for the neutral member within 
the same multiplet. ^+(4430) decays into 7r+-0' via S- 
wave. Moreover their scheme requires the ground state 
with J^*^ = 1^ around 3880 MeV which decays into 
TT^'ip and r]cP~^. 

With a QCD-string model, Gershtein, Likhoded and 
Pronko argued that both X(3872) and Z+(4430) are 
tetraquark states [6] . They speculate that the two quarks 
and two anti-quarks sit on the four corners of a square 
while any qq pair is a color-octet state. The decays of 
tetraquarks involve the reconnection of the color string. 

Cheung, Keung and Yuan discussed the bottom analog 
of Z+ (4430) assuming it is a tetraquark bound state [3] ■ 
According to their estimate, the doubly-charged Zt,c state 
lies around 7.6 GeV while the bottomonium analog Z^b 
of Z+(4433) is about 10.7 GeV. 

Qiao suggested Z+(4433) be the first radial excitation 
of Ac — S° bound state [8]. Within Qiao's scheme, all of 
the recently observed states r(4260), r(4361), Z+(4430) 
and y(4664) are accommodated in the extended heavy 
baryonium framework. 

Lee, Mihara, Navarra and Nielsen calculated the mass 
of Z+(4430) mz = (4.40 ± 0.10) GeV in the framework 
of QCD sum rules, assuming it is a 0~ molecular state 
of i:)*(2010)Z)i(2420) [i]. They predicted the analogous 
mesons Zg at mz^ = (4.70±0.06) GeV, which is above the 
D*Di threshold and Zbb around m^,, = (10.74 ± 0.12) 
GeV. 

Bugg proposed that Z+(4430) is a threshold cusp aris- 
ing from the deexcitation of the D*(2010)Di(2420) pair 
into lower mass D states [10|. The imaginary part Imf{s) 



of the elastic S-wave amplitude /(s) is a step function 
near the threshold. From the dispersion relation, the 
real part of the amplitude near the threshold looks like 
Ref{s) ~ /^^ ^ds' - In(s-so). Therefore |/(s)|2 con- 
tains a sharp cusp near threshold. 

However, none of the above schemes explains why 
Z+(4430) does NOT decay_into ir+ip. One notes that 
the momentum of D and D is small and close to each 
other in the rest frame of the parent B meson, especially 
when one (or two) of the D meson pair is an excited state. 
Therefore, there is plenty of time for the D meson pair 
to move together and re-scatter into 7r+-0'. However, the 
most puzzling issue of the re-scattering mechanism is the 
absence of any signal in the ■k'^ J/ip channel. One may 
wonder whether the mismatch of the Q-values of the ini- 
tial and final states plays an important role. If so, one 
should also expect a signal in the 7r"'"-0(3S') channel since 
there is nearly no mismatch of the Q-value now. Another 
potential scapegoat is the specific nodal structure in the 
wave functions of the final states. Detailed calculations 
along the above two directions are highly desirable to in- 
vestigate the origin of the non-observation of Z+(4430) 
in the tt+V' mode. 

In the heavy quark limit, the angular momentum of the 
light quark ji — I + Sq is a good quantum number where 
I is the orbital angular momentum and Sq is the light 
quark spin. For the P-wave heavy mesons, J; = | or ^ 
which correspond to the two doublets with J^ — (1+, 2+) 
and (0+, 1+) respectively. The ground states D and D* 
belong to the J^ — (0~, 1~) doublet with ji = i. Since 
the l"*" state in the (1+, 2+) doublet decays into D*Tr via 
D-wave, it's very narrow and denoted as I?i(2420) [ll|. 
The 1+ state in the (0+, l"*") doublet decays into D*7r via 
S-wave. Hence it's very broad and denoted as £''^(2430) 

Assuming Z+(4430) is a DiD* (or D[D*) S-wave reso- 
nance, Meng and Chao found that the open-charm decay 
mode D*D*7T is dominant and the re-scattering effects 
are significant in DiD* channel but not in D'lD* chan- 
nel since D'^ is very broad [ij]. For the J^ = 1^ can- 
didate, the ratio t{z+ -^ ij'Tr+)/r{Z+ -^ J/-07r+) may 
reach 5.3 with a special set of parameters, which partly 
accounts for why the Z+(4430) is difficult to be found in 
J/V'7r+. 

Despite so many theoretical speculations proposed 
above, a dynamical study of the Z+(4430) signal is still 
missing. In this work we will explore whether Z"'"(4430) 
could be a S-wave molecular state of D* and D[ (or -Di), 
which is loosely bound by the long-range pion exchange 
potential. We want to find out whether there exists the 
attractive force between D* and D[ (Di) in different 
channels. 

This paper is organized as follows. We discuss the 
possible quantum numbers of Z+(4430) and its possible 
partner states in Section [Hi We collect the effective La- 
grangians and various coupling constants in Section IlIII 
We derive the one-pion-exchange potential (OPEP) in 
section HVl Then we present our numerical result and a 



short discussion in Section FVl 



II. QUANTUM NUMBER OF Z+(4430) AND 
OTHER POSSIBLE STATES 



Natively, the smaller its angular momentum, the lower 
the mass of Z+(4430). Since Z+(4430) fies very close to 
the i:)*(2010)£'i(2420) threshold, we consider only the 
possibility of Z+(4430) being the loosely bound S-wave 
state of D* and D[ (or Di). Therefore, its possible an- 
gular momentum and parity are J^ = 0~, 1~ , 2~ . More- 
over, Z+(4430) was observed in the ^'tt"*" channel. So it 
is an isovector state with positive G-parity, i.e., J*^ = 1+. 

For a charged member within a molecular isovector 
multiplet, one can construct its flavor wave function with 
definite G parity in the following way. Suppose |^) is one 
component of its flavor wave function. Then the G = + 
state reads: 



\+) = ^(\A) + G\A) 



(1) 



where G = e^^^^C is the G-parity operator. Similarly, 
the G — — state reads: 



\-)^j={\A)-G\A)) 



(2) 



In the present case, the flavor wave function of 
Z+(4430) is 



\Z+) = ^[\A') + \B')] 



(3) 



with \A') = \D'I'D*+) and \B') = \D'+D*^), where D'^ 
and D'^ belong to the (0+,l+) doublet in the heavy 
quark effective field theory. Or 



\Z+) = —[\A) + \B) 



(4) 



with 1^) = |D0£I*+) and \B) = \D+D*^), where Dj and 
D-^ belong to the (l+,2+) doublet. 

The flavor wave functions of the partner states of Z+ 
can be derived in the following way. Charmed mesons 
belong to the fundamental representation of flavor SU{'d). 
Therefore, the system with a charmed meson and an anti- 
charmed meson belongs to 3x3 = 8-t-l. We list the wave 
functions of these hidden charm states, whose names are 
listed in Fig. [TJ Here we use the system of DJ-. and £*(«)! 
as an illustration. 

The flavor wave functions of these states are 
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der charge conjugate transformation, we have 



® 



FIG. 1: The multiplets composed with charmed mesons and 
anticharmed mesons. 



CJ2C-1 = -(c'^7^75d'^)(d^7^c'^) = -Ji. 

CJ4C-I = -{C^^,,J5U''){U^J^C^) = -J3 ■ 



Therefore, we get 



CJzoC 



'-[-J2+cJl-{-J4 + cJ3)]. 



(6) 
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d;,d:+)-c\d:-d:. 



We use Z° in the J = case as an example to iUustrate 
how to determine c. At the quark level, this molecular 
state may be written as 



Jz» = -^[Ji- c-h - (Js - CJ4)] 



(5) 



In other words, the C-parity oi Z^ is C = ± for c = ±1. 
Since Z+(4430) is a state with isospin 1 and G = +, one 
requires c — ~l for these states. The quantum numbers 
of Z° are I^{JP^) = l+(0, 1, 2) — . 

We want to emphasize that the presence of both the 
charm and anti-charm quark (and the light quark and 
anti-quark) in the expression of Ji ensures there is no 
arbitrary phase factor under charge conjugate transfor- 
mation. 

There is one intuitive and natural way to interpret 
the above results if we consider the flavor SU(4) sym- 
metry, which is of course broken badly in reality. If we 
naively assume the flavor SU(4) symmetry, then states 
within the same multiplet should carry the same coef- 
ficient under charge conjugate transformation. For ex- 
ample, C\D*-) = C{D*-)\D*+) where the coefficient 
C{D*~) takes the same value as either C{p°) or C{J/tp)- 
I.e., C{D*-) = -1. Similarly, C{D:[) = +1. In other 
words, the c — —I case leads to the Z° state with nega- 
tive C-parity. 

In addition, one obtains the flavor wave functions of 
these states with opposite G-parity if we take c = -f in 
the above equations. For example, we will also discuss 
whether Z+ could be a molecular state: 



|Z+) = ^[\D'lD*+)-\D*'^Dt 



If Z+ exists, this state may be discovered in either 
J/V'tt+tt'^ or iP'tt^tt'^ channel. 

If we replace one of the c (or c) by b (or 6), we can 
get molecular states such as (bq) — (cq). With the dual 
replacement c ^f h,c ^ b, we get the hidden bottom 
molecular states [bq] — (bq)- 



where 



Ji = (c"7M75d")(d^7^c^), 
J2 = (J''7M75c")(c^7^d'), 
h = (c'^7m75u'^)(S^7^c^), 
J4 = (w'^7p75c'^)(cVw')- 

In the above equation, a and e are the color indices. Un- 



III. EFFECTIVE LAGRANGIANS AND 
COUPLING CONSTANTS 



We collect the effective chiral Lagrangian used in the 
derivation of the OPEP in this section. In the chiral and 
heavy quark dual limits, the Lagrangian relevant to our 



calculation reads [Ij, [lj| 

C = igTr[Hb4bal5Ha] + ig'TT[Sb4baj5Sa] 

+ig"Tv[Tf,i,4bal5T^] 

+ [ihTi[Sb4baj5Ha] + h.c] 






.{,^T,[T,^(^]pA,)ba75Ha] + h.c.}, (7) 



where 



Ha = ^[Pr-Pa75], 
Sa = ^[P['al,l5-PLh 



rpfl 






(8) 
(9) 

(10) 



and the axial vector field A^j^ is defined as 

Kb = \{^^d'^ - id^'i^U = j^d^M + • • • 

with ^ = exp(?;Al//^), U = 132 MeV and 

M=\ .- --; + ^ i^o . (11) 

K- Ro -^ ; 

After expanding Eq. ([7]) to the leading order of the 
pion field, we further obtain 

(12) 



£ 



.'OD'o^oe'^^"'D[l{d,7T°)id^D%) + h.c 



d'°d["t!o — 9d[°d[°tt'>'' ■ -^iqV^m" yv^i^^i/s; 



(13) 

>Cl50^o,o = too^o,oe"''^''(7?°J(-9p^")(a.50^) + /i.C., 

(14) 

+D;{d'^7T°){d^D[j + {d^D:){d,7T°)D["] 
+h.c. (15) 



^D*Di7rO = gD'Di-no[D*pP'lvg^''{d'^d\'K°) 



doublet. The coupling constants gD'+D'+Tr^-, 9d'°d''>t:°^ 
9d°d°7,o and g^.^O^o are 



_V2g' 



V2g 
5.9" 



_ _ i\/2/i 



^mD'TTlDi 



gn^ono.o--g^.,^.^o--^^j-^^ 



{hi + h2). 



-DlDi^g^'^id^dxTT^) + 2D;Di^g^\d''dx7r°) 
+ ^ {d^D*n{dadxTr''){d"Di,)], (16) 



mD'TUDi 



where DJ^ denotes the P-wave axial-vector state in the 
(0+, 1+) doublet while Di is the 1+ state in the (1+, 2+) 



The coupling constant g was studied in many theoret- 
ical approaches such as QCD sum rules [H, [la, lOj [3 
and quark model [l3|. In this work, we use the value 
g = 0.59 ± 0.07 ± 0.01 extracted by fitting the experi- 
mental width of D* jJJi] . Falk and Luke obtained an ap- 
proximate relation \g'\ = \g\/3 and \g"\ = \g\ in quark 
model [l3|. However the phase between g' and g" is 
not fixed. With the available experimental information, 
Casalbuoni and collaborators extracted h = —0.56 ±0.28 
and h' = {hi + h2)/A^ = 0.55 GeV"! [H. If we re- 
place the meson field in the heavy quark limit in the 
above equations by the fields in full QCD and scale the 
coupling constants by a factor ^vhb' etc, we get the ef- 
fective Lagrangian in full QCD, which is used below in 
the derivation of the potential. 



IV. DERIVATION OF THE ONE PION 
EXCHANGE POTENTIAL 

Study of the possible molecular states, especially the 
system of a pair of heavy mesons, started more than 
three decades ago. The presence of the heavy quarks 
lowers the kinetic energy while the interaction between 
two light quarks could still provide strong attraction. 
Okun and Voloshin proposed possibilities of the molecu- 
lar states involving charmed mesons [2(|. Rujula, Geogi 
and Glashow once suggested -0(4040) as a D*D* molecu- 
lar state |2l|. Tornqvist studied possible deuteronlike 
two- meson bound states such as DD* and D*D* us- 
ing a quark-pion interaction model |22l | . Dubynskiy and 
Voloshin proposed that there exists a possible new reso- 
nance at the D*L»* threshold [2l,[2i|. 

Several groups suggested X(3872) could be a good 
molecular candidate [2^, l2y, [23, [2^, [23| . However, Suzuki 
ajgued that X(3872) is not a molecule state of D°D*° + 
D°D*° ^. Instead, X(3872) may have_a dominant 
cc component with some admixture of D'^D*^ + D^D*^ 
30, 31, 32]. 

In this work we will explore whether Z+ (4430) could be 
a S-wave molecular state of D* and D'l (or Di), which is 
loosely bound by the long-range pion exchange potential. 
We first derive the scattering matrix elements between 
the pair of D mesons as shown in Fig. [51 With the Breit 



approximation, we can get the one-pion exchange poten- 
tials. Since the flavor wave function of Z"*" contains two 
components, we have to consider both the direct scatter- 
ing diagram Fig. [2] (a) and the crossed diagram Fig. [2] 
(b). Note only the crossed diagram contributes to OPEP 
in the case of X (3872). 

Recall the flavor wave function of Z^ (4430) reads 



|Z+) = i=[|A') + |S'> 



with \A!) = \DfD*+) and \B') = \D'+D*°) or 



(17) 



(18) 



with \A) = \DlD*+) and \B) = \D+D*"). The mass of 
Z+(4430) is expressed as 



Mz+ = mo- + m ,>) +T + E + 5, 

^1 



(19) 



where T is the kinetic energy in the center of 
mass frame, E = {D*+Dp"\ni\D*+Dp°) and 6 = 

{D*+Dp°\n2\D'+D*°). Hi and 0.2 correspond to the 
interaction in Fig. [2] (a) and (b) respectively. For the 
possible Z^ state with negative G-parity, we can get its 
mass through the replacement +S -^ —S in Eq. (|19p . 
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The single pion exchange in the direct scatter- 
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A. Scattering amplitudes 

We collect the scattering amplitudes in the 
different channels below. For the process 



D*+(pi,ei).D;o(p2,e2) ^ i5*+(p3, e3)5'i°(P4, £4), the 
amplitude is 



s'o 



M[D*+{p,,e,)Dl\p2,e2) -^ D*+{p^,e^)Dl\pi,ei)] 



^'^/ ,_^ ^ ^"''^'^^"'^'"'"'gMVPi^Pi^' 



/2 gr2 _ ^2 



(20) 



For D*+{pi,ei)D\{p2,e2) -> D*+{p3,e^)D\{pi,ei), one 
gets 

iM[D*+{p^,ei)Dl{p2,e2) ^ D*+(p^,es)Dl{pi,ei)] 
bigg" 1 



3/2 92 „ ^2 

^('^1q'=2q')(£3^^4^')- 



e''^^''e'''P'^''''q^q^,p^,Pi,, 



(21) 



For i?*+(pi,ei)5;o(p2,e2) ^ i?'+(p3, 63)0*0(^4, £4), the 
amplitude is 

iM[D*+[p^,e^)D[\p2,e2) ^ i?;+(p3, e3)5*°(P4, £4)] 



2i/i2 1 



[-P?/ - g>!^ + (Pi • 9)5"^] 



a/3l 



/2 g2_^2i 

xhp?'/' - 9" pf + (P4 • q)9^'^']{4h4i')i4i4}')- 

(22) 

The amplitude of the process £'*+(pi, ei)Z5]'(p2, £2) ^ 
I?+(p3,e3)5*°(p4,e4)is 

iAl[i?*+(pi,ei)-D°(p2,e2) ^ i?^(p3,e3)^*°(p4,e4)] 
_ —imD*mDi{hi + h2)^ 1 



3 



/3iy 






3/.2A2 

-(pi ■q){q-p3] 



q2 _ ^2 



3g'^g'3 - /''g^ 



3//'-/''^'g2 






(23) 



Here the polarization vector is defined as e*^ = 
^(0, ±1, i, 0) and eO = (0,0, 0,-1). 

B. The one-pion-exchange potential 

We impose the constraint on the scattering amplitudes 
that initial states and final states should have the same 
angular momentum. The molecular state \J,Jz) com- 
posed of the 1~ and 1"*" charm meson pair can be con- 
structed as 

|J,J,)- ^(l,Ai;l,A2|J,J,)bi,ei;p2,e2) (24) 

Ai,A2 

where (1, Ai; 1, A2I J, Jz) is the Clebsch-Gordan coeffi- 
cient. Combining the equation with the scattering am- 
plitudes, one gets the matrix element iM.{J, Jz). 



With the Breit approximation, the interaction poten- 
tial in the momentum space is related to iM{J, Jz) 



V{q) = 



1 



n, 2m, H/ 2to/ 



-.M{J,Jz) (25) 



where mi and mf denote the masses of the initial and 
final states respectively. We collect the expressions of the 
potential in Tables [J and [ill We have explicitly shown 
the resulting potentials are the same for the different Jz 
component. 



TABLE III: The one pion exchange potential in the coordinate 
space with A' = Of D*+ and B' = D'+ D*'^ . 



State 



A'{B') -^ A'TW) 



99 Ufr) - HiiLg- 



fj[^(r) - 



SS-[5(r) 



6/ 



A' {B') -* B' (A^ 



h-'(q")'-' cos(Mr) 



Stt/I 



h-'(q")-' cos(Mr) 



TABLE I: The one-pion-exchange potential between D* {D*) 
and Di (Di). Here the expressions are for the J^ = compo- 
nents. A' = D[°D*+ and B' = D'+D'°. 



State 


A'{B') ^A'{B') 


A'{B') ^ B'{A') 


0" 


' '2 

gg q 

3fi q^+mj 


2fi (q"y'-^'-m-i 


1" 


2f4 q^+mi 


h^ (g'Y 
2fi (q'O^-q^-™^ 


2" 


gg' Zq^'-agf 





TABLE II: The one-pion-exchange potential between D* {D* 
and Di (Di). In this table, A = I)?D*+ and B = D+D*°. 



State 


A{B) ^ A{B) 


A(B) -» B(4) 


0" 


ngg" q^ 

ISfi q^+mj 


(if'i (9'')^-q^-mi; 


1" 


599" 9j 
12 f^ q-'+m-i 


h'2 2{q^)^-39^q^ 
12/^ (q.l)2_q2_™j 


2" 


599" 2q^-3<J^ 
36/^ q^+mj 


h'2 (q-^)^_8g^q-'+99; 

8/^; (9")^-q^-»"i 



1. OPEP from the direct scattering diagram 

From Tables H] and HIl it's very interesting to note 
that the OPEP from the direct scattering diagram in the 
J^ ~ 0" and 1~ channel always has the same sign while 
it's opposite to that in the J^ = 2~ channel. In other 
words, there must be attraction in one of three channels 
no matter what sign gg' takes. 

After making the Fourier transformation, we get the 
one-pion-exchange potential in the configuration space. 
The final potentials are obtained with the following re- 



placements: 



73, 



75: 



9 ? 

X y 



r /lb etc 



since we consider only the S-wave system. Alternatively, 
one may average the potential in the momentum space 
first. 

In the derivation of the OPEP, we let g" w as usually 
done for the direct scattering diagram in Fig. [21(a). The 
resulting potential from the direct scattering diagram is 
the familiar Yukawa potential plus a 5 function, similar 
to OPEP in the nucleon-nucleon potential. 



TABLE IV: The one pion exchange potential in the coordinate space with A = D^D* and B — DfD* . 



State 



AiB) ^ AiB) 



18£ 



■[S{r) 



' ^99" \x(j.] _ Hl^g- 



599 Ufrl - ^^e" 



AiB) ^ BiA) 



6£ 



-[\/'5{r)~^,^S{r) 



h'2 

US 



[V^5(r) - M'5(r) . 



fj. cos£ir 1 
47r r J 



60/i 



.[\/^S{r)-i,'S{r) 



r 



2. OPEP from the crossed diagram 

First we note from Table [TTl that the OPEP from the 
crossed diagram contains terms such as q^ in the numer- 
ator, which reflects the fact that Di decays into D*tt 
via D-wave. But the approximation q° « is not rea- 



sonable for the crossed diagram in Fig. [2l (b) because 
q° Ri Mdi - Md* « 410 MeV « 3m^. In the following, 
we take J^ — 0^ as an example to illustrate how to deal 
with the crossed diagram. Since the D[ meson decays 
into Z?*7r, the crossed diagram contains a small imagi- 
nary part, which is roughly of the order T^)' . However 
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FIG. 3: (a) The solid, dotted and dashed hnes correspond 
to the potentials from the direct scattering diagram A(B) -^ 
A(B) in the J^ = 0~,1~,2~ channels respectively, where 
g — 0.59, g' = g/3 and g" = g. (b) The potential from the 
crossed diagram A{B) -^ B{A). 



(b) 

FIG. 4: (a) The solid, dotted and dashed lines correspond to 
the potentials from the direct scattering diagram A'(B') -^ 
A{B') in the J^ — 0~,1~,2~ channels respectively, where 
g = 0.59, g' = g/3 and g" — g. (b) The potential from the 
crossed diagram A'{B') -^ B'{A'). 



only the real part of the scattering amplitude contributes 
to the potential. Hence the principal integration is always 
assumed throughout the Fourier transformation for the 
crossed diagram. 



V{v) 



(g' 



,0^2 






2/2(g0)2-q2-m2+^e 
(9°)' „,o.rl dd 



o«qr 



dq 



o«qr 



/I- 



(2^) 



h {q^y cos(^r) 

where ^^ = (q^)^ - m\ and g° sa Md^ ~ Md*- The 
resulting potential is of very long range and oscillating. 
We list the expressions of the potential in the coordinate 
space in Tables Hill and IIVI 

For the D*D^ system, it's clear from Table Hill that the 
OPEP from the crossed diagram has the same overall sign 
in all three channels. Moreover, the sign is positive. Be- 
cause of the factor Qq ~ 9m^ , the OPEP from the crossed 
diagram is numerically much larger than that from the 



direct scattering diagram. For the D*Di case, we note 
from Table ITVl that the overall sign in the 1~ channel is 
different from that in both 0~ and 2~ channels. 

With the couphng constants g — 0.59, g' = g/3 and 
g" = g, the variation of the potential listed in Tables Hill 
and IIVI with r (in unit of GeV^^) is shown in Figs. [3] and 

SI 

We also note that the OPEP does not depend on the 
mass of the charmed mesons. In other words, the same 
potential may be used in the discussion of the hidden 
bottom molecular states (bq) — (bq). 



RESULTS AND DISCUSSIONS 



Besides those coupling constants in Section HUl we also 
need the following parameters in our numerical analysis: 
mo' = 2007 MeV, mD[ = 2430 MeV, mo, = 2420 MeV, 
TOB* = 5325 MeV, ms; = 5732 MeV, j^ = 132 MeV, 
m^ = 135 MeV [Hi; rriB, = 5725 MeV [H. 

With the potentials derived above, we use the vari- 
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ational method to investigate whether there exists a 
loosely bound state. Our criteria of the formation of 
a possible loosely bound molecular state is (1) the ra- 
dial wave function extend to 1 fm or beyond and (2) the 
minimum energy of the system is negative. Our trial 
wave functions include (a) ipi'>') = (1 + ar)e~°"'; (b) 
^(r) = (1 + ar2)e-'''''; (c) VW = ^^(1 + ar)e-^^ 

Unfortunately a solution satisfy the above criteria does 
not exist for the system oi D^~D* or Di^D* in all J^ = 
0~, 1~, 2~ channels with the realistic coupling constants 
g — 0.59, g' = g/3 and g" — g deduced from the width of 
D*,Di and D[. Such a solution also does not exist if we 
switch the sign of gg' or enlarge the absolute value of gg' 
by a factor 3. The same conclusion holds for the system 
oi b' - B*, B - B* and Z+ with negative G-parity. 

In short summary, we have performed a dynamical 
study of the Z^ (4430) signal to see whether it is a loosely 
bound molecular state of Di — D* or D[ — D* . We find 
that the interaction from the one pion exchange potential 
alone is not strong enough to bind the pair of charmed 
mesons with realistic coupling constants. Other dynam- 
ics is necessary if Z+(4430) is further established as a 
molecular state by the future experiments. 

It's interesting to note that the one pion exchange 
potential alone does not bind the deuteron in nuclear 
physics either. In fact, the strong attractive force in 
the intermediate range is introduced in order to bind 
the deuteron, which is sometimes modeled by the sigma 



meson exchange. One may wonder whether the simi- 
lar mechanism plays a role in the case of Z+(4430) and 
X(3872). Further work along this direction is in progress. 
If the Z+(4430) is really a J^ = 1" molecular 
state, the quantum number of its neutral partner state 
Z'''(4430) is J^*^ = 1^^. Such a state can be searched for 
in the e^e~ annihilation processes. Babar and Belle col- 
laborations have observed several new charmonium (or 
charmonium-like) states with J^'^ = 1 around this 
mass range including y(4260), y(4320), and r(4664) 
with the initial state radiation (ISR) technique, although 
these states do not appear as a peak in the R distribu- 
tion. We strongly urge Belle and Babar collaborations to 
search for the Z°(4430) state in the ■k'^iJ}' channel using 
the ISR technique. The absence of a signal will be an 
indication that the J^ of Z+(4430) is not l~ . 



Acknowledgments 

We would like to thank Professor K.T. Chao, Professor 
Z.Y. Zhang, and Dr C. Meng for useful discussions. This 
project was supported by the National Natural Science 
Foundation of China under Grants 10421503, 10625521, 
10675008, 10705001, 10775146 and the China Postdoc- 
toral Science foundation (20060400376). 



[1] Belle Collaboration, K. Abe et al., arXiv:0708.1790] [hep- 

ex]. [is; 

[2] Barbar Collaboration, B. Aubert et al., Phys. Rev. D 71, [19 

071103 (2005). 
[3] Barbar Collaboration, B. Aubert et al., Phys. Rev. D 73, [20 

011101 (2006)_ 

[4] J. L. Rosner, [MXiv:0 708.3496 [hep-ph]. [21 

[5] L. Maiania, A.D. Polosab and V. Riquerb, 

arXiv:0708.3997 [hep-ph]. [22 

[6] S.S. Gershtein , A.K. Likhoded and G.P. Pronko, 

larXiv:0709.2"058 [hep-ph]. [23 

[7] K. C heunga, W. Y. Keung and T. C. Yuan, [24 

larXiv :0709. 1312 1 [hep-ph]. 

[8] C.F. Qiao, arXiv:0709.4066 [hep-ph]. [25 

[9] S. H. Lee, A. Mihara, F.S. Navarra and M. Nielsen, [26 

[arXlV :0710.1029 [hep-ph]. [27 

[10] D. V. Bugg, .arXiv:0709.1254 [hep-ph]. [28; 

[11] W. M. Yao et al.. Particle Data Group, J. Phys. G 33, 1 

(2006). [29 

[12] C. Meng and K. T. Chao. larXiv:0708.4222 l [hep-ph]. [30 

[13] A. F. Falk and M. Luke, Phys. Lett. B 292, 119 (1992). [31 

[14] R. Casalbuoni, A. Deandrea, N. Di Bartolomeo, F. Fer- 

uglio, R. Gatto and G. Nardulh, Phys. Rept. 281, 145 [32 

(1997). 

[15] V.M. Belyaev, V. M. Braun, A. Khodjamirian and R. [33 

Ruckl, Phys. Rev. D 51, 6177 (1995). 
[16] F. S. Navarra, Marina Nielsen and M.E. Bracco, Phys. 

Rev. D 65, 037502 (2002). 
[17] F. S. Navarra, M. Nielsen, M. E. Bracco, M. Chiapparini 



and C.L. Schat, Phys. Lett. B 489, 319 (2000). 

Y. B. Dai and S. L. Zhu, Eur. Phys. J. C 6, 307 (1999). 

C. Isola, M. Ladisa, G. Nardulh and P. Santorelli, Phys. 

Rev. D 68, 114001 (2003). 

M. B. Voloshin and L. B. Okun, JETP Lett. 23, 333 

(1976). 

A. D. Rujula, H. Georgi and S. L. Glashow, Phys. Rev. 

Lett. 38, 317 (1977). 

N.A. Tornqvist, Nuovo Cim. A 107, 2471-2476 (1994); 

Z.Phys. C 61, 525-537 (1 994). 

M. B. Voloshin, |arXiv:hep- ph/0602233 

S. Dubynskiy and M. B. Voloshin, Mod. Phys. Lett. A 

21, 2779 (2006). 

F. E. Close, P. R. Page, Phys. Lett. B 578, 119 (2004). 

M. B. Voloshin, Phys. Lett. B 579, 316 (2004). 

C. Y. Wong, Phys. Rev. C 69, 055202 (2004). 

E. S. Swanson, Phys. Lett. B 588, 189 (2004); ibid B 598, 

197 (2004). 

N. A. Tornqvist, Phys. Lett. B 590, 209 (2004). 

M. Suzuki, Phys. Rev. D 72, 114013 (2005). 

S. L. Zhu, hep-ph/0703225 , to appear in Int. J. Mod. 

Phys. E; arX iv:0707. 4586 [hep-ph], 

C. Meng, Y. J. Gao and K. T. Chao, arXiv: 

hep-ph/0506222l, 

CDF Collaboration, Mass and width measurement of 

orbitally excited (L = 1) B"° mesons, CDF Note 8945, 

www-cdf.fnal.gov/physics/new/bottom/070726.blessed- 

bss/; CDF Collaboration, Andreas Gessler, the EPS 

talk given on HEP 2007 in Manchester, .arXiv:0709.3148. 



[hep- 



